Scattering and binding of different atomic species in a one-dimensional optical lattice 
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The theory of scattering of atom pairs in a periodic potential is presented for the case of different 
atoms. When the scattering dynamics is restricted to the lowest Bloch band of the periodic potential, 
a separation in relative and average discrete coordinates applies and makes the problem analytically 
tractable. We present a number of results and features, which differ from the case of identical atoms. 

PACS numbers: 03.75.Lm, 34.10.+X, 63.20.Pw, 71. 23. An 



The combination of external confinement and tuning 
of atomic levels and molecular potential curves by exter- 
nal fields has led to a rich variety of experiments with 
cold atoms, aiming at the study of formation of ultra- 
cold molecules, mean field dynamics, generation of com- 
plex many-body states, and production of quantum cor- 
related sources of atoms. While initial studies dealt with 
only a single atomic species, the achievement of sym- 
pathetic cooling, mixing and observation of spatial sep- 
aration of bosonic and fermionic mixtures, and forma- 
tion of heteronuclear molecules with permanent dipole 
moments have spurred activities on mixtures of atomic 
species fliaiSllllIS]. 

In this paper, we consider the physics of a pair of atoms 
labeled A and B in an optical lattice consisting of three 
standing wave laser beams intersecting at right angles 
with wavelength A l • The atoms experience light induced 
energy shifts, which result in a cubic periodic potential 
'i^iat,/3(x) = X]-_;l,2,3'K^3^^^^(^^^/^)' where a = Xl/2 is 
the lattice constant and is the lattice strength for 
the atomic species P = B along direction i. The lat- 
tice strengths depend on the polarizability of the atoms 
as well as the amplitude and detuning of the electromag- 
netic field. In our numerical examples we apply the phys- 
ical parameters relevant to ^^K and ^''Rb atoms, studied 
in Ref. [1 . We take A to be ^'^Rb and B to be "^^K with 
Vl = 3^R,Rb, Vl = Vl= 40ER,Rb and = 0.86^1, 
where £^R,Rb = /i^/2mRbAL is the recoil energy for ru- 
bidium. 

The motion of a single atom is governed by the Hamil- 
tonian Hq = — (fi^/2m/3)V^ + Mat,/3(x), which separates 
into three independent one-dimensional equations, each 
solved by Bloch wave functions where n is the 

band index and q G [— 7r/a, 7r/a] is the quasi- momentum. 
Within each band the Bloch waves can be transformed 
into the localized Wannier basis functions w^^.{x) cen- 
tered around the potential minima zj = ja [6 . We are 
interested in the quasi-one-dimensional regime, which is 
obtained when the lattice potential in one (longitudinal) 
direction is much weaker than in the two other (trans- 
verse) directions <^ Vp = V^. In this situation 
the transverse motion will be confined to the ground 
state of an effective harmonic oscillator with frequency 



ujf3 = (7r/a)y 2y^/m/3 and eigenstate of motion w^{x) = 

{nif^cOf^/Trh)^^'^ ex.-p{— 171(300 f^x'^ /2h) whenever the longitu- 
dinal motional energies are significantly below huof^. The 
motion is then effectively one-dimensional along the di- 
rection Xs. 

The dynamics of a single atom in a one-dimensional 
lattice is described by the dynamical tunneling ampli- 
tudes Jp = {uJnzj\^o\^nzj+i)- following wc ouly 
include nearest-neighbor tunneling, and we consider mo- 
tion restricted to the lowest Bloch band (n = 1), but our 
analysis and results may be generalized to higher bands 
and beyond nearest-neighbor tunneling [7 . 

We characterize the two-body system by 
a wavefunction, which is expanded on Wan- 
nier product wavefunctions ^(^3,a,^3,b) = 
T^ZA^ZB ^i^A^^B)w^z^{xs,A)w^z^{xs,B) and solve 
for the amplitude '^(^a, ^b) of finding atom A and atom 
B in the Wannier functions centered at the discrete sites 
za and zb^ respectively. In this basis the two-body 
non-interacting Hamiltonian becomes 

Ho = -Ja{A,^ + 2) - Jb(A,, + 2), (1) 

where Azf{z) = f {z -\- a) -\- f {z — a) — 2f (z) is the discrete 
Laplacian. The eigenstates of Hq are product states 

iIj{za, zb) = exp{iqAZA) exp{iqBZB) (2) 
with energy 

<Qa.Qb) = EA{qA)^EB{qB). (3) 

where Ep{qp) = —2Jpcos{qpa) is the single particle en- 
ergy dispersion. For convenience we put the zero of en- 
ergy in the middle of the first Bloch band. 

The eigenstate in Eq. ([2| is a product of Bloch waves 
with quasi-momenta qA and qB- We now change to col- 
lective Z = {za-\- zb)/2 and relative z = za — zb coordi- 
nates, and we introduce the collective K = qA -\- qB and 
relative q = {qA — qB)/2 quasi-momenta. In the special 
case of equal masses the collective coordinate is identi- 
cal to the center-of-mass coordinate. Due to the discrete 
nature of the problem the Hamiltonian separates into a 
collective and a relative coordinate part via the Ansatz 
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iP{za,zb) = e^^^ iIjk{z). By applying to the product 
Ansatz we get Hoe'^^i/jxlz) = e'^^ H^iljxiz), where 
Hj^ is the action of the Hamiltonian Hq on the relative 
coordinate part ipx given by 

HUKiz) = -Ja (e*^"/VK(z + a) + e-'^'^^^iPKiz - a)) 

-Jb (e-^^°/2^K(^ + a) + e'^«/Vx(^ - a)) • 

(4) 

Note that the relative motion Hamiltonian acts on 
components of the joint system, where K has a specified 
value. 

Unlike the case of identical atoms, where Ja = Jb^ 
is not invariant under complex conjugation, and hence is 
not time-reversal invariant, for 7^ 0. The collective 
quasi-momentum is here playing a role similar to that of 
a classical magnetic field on the motion of electrons in 
an atom or a solid. As a consequence of the breaking of 
the time-reversal symmetry, we cannot expect the bound 
state wavefunctions to be real, but in further analogy 
with the magnetic interactions we note that time rever- 
sal of the full two-body dynamics is accompanied by a 
change of sign of K (of the magnetic field), and hence 
= {H^j^y. If tp^K is an eigenstate of i^^^ then 
the complex conjugate, and hence time-reversed, wave 
function V^* ^ is an eigenfunction of H^^ with the same 
eigenvalue. 

Introducing the average and half-difference tunneling 
amplitudes, Js,a = {Ja^Jb)/"^^ with values Jj: = Ja = 
Jb and Ja = in the case of identical particles, the 
energies (|3| can be written in the convenient form 



exiq) = Ek cos [(Tk + q)a] , 
which is parametrized by the collective energy 



(5) 



Ek = -'^\JJl cos2(ifa/2) + J2 sm^{Ka/2) (6) 



and by the quasi-momentum shift 

_i 4Ja sm{Ka/2) 



Tfca = sin 



\E 



(7) 



where T^a is taken to be in the interval (— 7r/2,7r/2] for 
Ka G (— 7r,7r]. The collective energy determines 
the extremal values of the band of energies obtained by 
varying g, and therefore the total width of the band is 
2\Ek\' The continuum band is shown as a the colored 
region in Fig. [ija). For the Rb-K system and our choice 
of optical lattice parameters the tunneling amplitudes 
have the values Jj: = O.SkHz and Ja = O.SJ^]. 

There are two significant differences from the identical 
particle case: {i) the width of the continuum band does 
not vanish for = 7r/a, but maintains a finite width of 
8| Ja|. For atom pairs with K = the width of the band 
is 8| J5]|. (ii) The minimum and maximum of the contin- 
uum are obtained for relative quasi-momenta q = —Tk 
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FIG. 1: (a) Energy band spanned by exiq) for the ratio Ja = 
0.5Je. The single atom spectra E(3{q(3) in (b) and the two- 
body spectrum exiq) = EA{K/2-\-q) -\- Eb{K/2 — q) in (c) are 
shown for K — 0.87r/a. The corresponding energy minimum 
{K, q) = (0.87r/a, —Tk) is marked by dots in all three figures 
and the vertical lines in (b) indicate from left to right qA, K/2 
and qB- The lower and upper curve in (a) shows the bound 
state energy for U = — 2.7Js and U = 2.7Js, respectively, 
while the shading shows the reflection coeflficient for \U\ — 
2.7 Jy: [see Eq. ( |13[ )]. The signiflcance of the diamonds in (c) 
is explained after Eq. (12). 



and q = ±7r/a — T^, respectively, and not, as for identical 
particles, at relative quasi-momenta and ±7r/a. 

The finite quasi-momentum difference at the band 
edges, qA — qB = —2Tk (mod 27r/a) , appears because 
the energy dispersions of the atoms have different am- 
plitudes, 2 J/3. This point is illustrated in Fig. [ija)- 
(c), where the lowest energy state within the band for 
Kajn = 0.8 is illustrated by dots. The quasi-momenta 
of the individual atoms are displaced by ±Tx from i^r/2. 

The value of Tk is depicted in Fig. |2j Tk is an odd 
function of both K and Ja- In Fig. |2|a) the lower curve 
Ja/Jt. = corresponds to J a = Jbi the identical par- 
ticle case, where Tk = for all K because both atoms 
have the same energy dispersion, and the upper curve 
Ja/Jt. = 1 corresponds to Jb = 0, where only atom A 
is allowed to move, and therefore qA = qB = and 
Tk = K/2. From Fig. |2|b) we note that Tq is always 
zero, and whenever J a Jb we have T^/a = 7r/2a due 
to the symmetry of the single particle dispersions. 

So far, our analysis has only accounted for the sepa- 
rable state ([2| of non-interacting atoms in a collective 
set of coordinates [8]. In the presence of an atomic in- 
teraction depending only on the relative coordinate the 
separation in z and Z is maintained, and we now turn to 
the description of the dynamics and bound states of the 
interacting system. In the Wannier basis we assume the 
interaction potential, [/, to be on-site, i.e., of the form 
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FIG. 2: Shift of the relative momentum, Tk- (a) From below 
Ja/Js = 0, 0.1, 0.25, 0.5, 1. (b) From below Ka/n = 
0, 0.5, 0.8, 0.95, 1. 



U5z^o^ where is a function of the three-dimensional 
background scattering length and the optical lattice pa- 
rameters. The bound and scattering states are analyzed 
using the relative motion Green's function for the non- 
interacting particles, G%{E) = [E — ^^o]~^- Due to the 
assumption of on-site interactions, we only need to spec- 
ify G%{E,z) = {z\G%{E)\0) given by the Fourier trans- 
form 



ae" 



dq 

^/^2^E- Ek cos{{Tk + q)a) + ir] 



(8) 



of the relative quasi-momentum Green's function, 
Gk{E; q') = d{q — q')/{E — ex(g)). Here is a positive 
infinitesimal added to enforce outgoing boundary condi- 
tions. By a change of coordinate q ^ q — Tk the integral 
is identical, except for a front factor exp(— zT^^;), to the 
one studied previously for identical particles [9]. Hence, 
for energies inside the continuum band we have 



G%{E,z) = 



(9) 



where pa = cos ^{E/Ek)^ and outside the continuum 
we find 

G'k{E,z) = sgn{E f ^^!l^ [-sgn(E)]-/° (10) 



with tza = cosh~^ \E/Ek\' The Green's function is sim- 
ilar to the identical particle case, except for a complex 
factor e~*^^^, as already mentioned, and a modified ex- 
pression for Ek 

We use the Dyson equation, G^{E) = G^i^) + 
Gk{E)UG^{E)^ to find the interacting Green's function 



1-UG%{E,0)' 



(11) 



and the scattering wavefunction then follows from the 
Lippmann-Schwinger equation 

tljK{E,z) = e'^P-^^^' ^UG^{E,0)e-'^^'e'P\'\. (12) 



Note that the relative quasi-momenta of the incoming 
q = p — Tk and reflected q' = —p — Tk waves are not 
related in the usual way, since q' ^ — unless K = {). 
This is explained by the degeneracy of q and q' due to the 
symmetry of the eK{q) around the minimum energy state 
q = — Tx, as indicated by the diamonds in Fig. [ijc). 

From the scattered wave in Eq. ( p!2| ) we can identify the 
scattering amplitude f{E^K) = UGk{E^O) and thereby 
the reflection coefficient 



R{E,K) = \f{E,K) 



El-E-^ + m' 



(13) 



which is indicated by the shading in Fig. [TJa). It reaches 
unity at the continuum boundaries, where the density 
of state diverges, and its minimum at the center of the 
energy band. The transmission coefficient T{E^ K) = 
1 — R{E^ K) might be probed in scattering experiments 
or measured spectroscopically by radiative coupling of a 
bound molecular state to the continuum [9 . In the latter 
case a deep-lying, tightly confined molecular state \i) is 
coupled to the continuum by some transition operator T, 
and the transition probability is then given by 

\{^K{E)\f\i)\^ a \^k{E,z = {))\^ = T{E,K). (14) 

The system supports bound states of the atoms, found 
as the poles of the scattering amplitude, f{E^K)^ with 
energy E^ = sgn{U)y^E'^ + /7^. As for identical parti- 
cles, we find a bound state below the continuum in the 
case of attractive interactions and a repulsively bound 
state [10 above the continuum in the case of U > 0. 
The corresponding bound state wavefunction of relative 
motion is given by G^{E^^ z)^ Eq. (10), up to a normal- 
ization factor. This function is exponentially decaying 
with \z\. 

The bound state wavefunctions are complex with 
a spatial phase variation of the relative motion, c.f 



Eq. (10). This is in agreement with our earlier discussion 
of the lack of time-reversal invariance of the relative mo- 
tion Hamiltonian restricted to fixed values of K. While 
such a complex phase variation may give the impression 
that one particle is passing by the other one over and 
over again within the exponential envelope of their rela- 
tive motional state, the two atoms are actually moving 
with the same group velocity. To see this, we recall that 
the bound states above (below) the continuum have pre- 
dominantly the quasi-momentum components of the non- 
interacting continuum states close to the upper (lower) 
band edge. This implies that the quasi-momentum dis- 
tribution of the attractively and repulsively bound states 
are peaked around q = —Tk and q = ±7r/a — T^, respec- 
tively, which can also be recognized directly from Eq. 
([To]). For these eigenstates of Hq a simple calculation 



deKjq) _ dEAdqA ^ d^;^ dqs 
dq dqA dq dqs dq 



h{vA - vb) (15) 
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FIG. 3: Quasi-momentum probability distribution in units of 
(a/27r)^ for the species A and S in a bound state with (a) 
U = -2.7Js and (b) U = +2.7Je. 



relates the group velocities vp of species P to the deriva- 
tive of the energy dispersion. At the band edges exiq) 
reaches its extrema, i.e., deK{Q)dq = 0, exactly when 
q = —Tk or q = ±7r/a — T^, and the group veloc- 
ities agree for these relative quasi-momentum compo- 
nents. Furthermore, since exiq) is an even function 
around q = — T^, the derivative and hence the differ- 
ence in group velocity va — v^j is odd. If we return to 

we therefore observe 



(12), 



the scattering state in Eq. 
that the difference in group velocity always changes sign 
when the wave is reflected. 

In Fig. |3] we show the quasi-momentum distribu- 
tion of the two components A and B for both a repul- 
sively (a) and attractively (b) bound atom pair [11 . We 
have assumed a Gaussian distribution in K with center 
K = O.Sn/a and standard deviation ax = O.lir/a. Clear 
peaks appear dit qA -\- qs = K and qA — qs = — 2T^ as 
marked by crosses. Here = 0.327r/a, corresponding 
to (qA^qs) = (0.087r/a, 0.727r/a) in the attractive case 
(a), while (qA^qs) = (— 0.927r/a, — 0.287r/a) in the case 
of repulsion (b). 

We suggest measuring the wave number shift Tk by 
rapidly turning off the lattice potential, because an adi- 
abatic ramp down would alter Ja and thereby Tk- The 
shift of the main peak in the momentum spectrum of each 
individual species is then given by K/2 K/2 ± Tk for 
[/ < and K/2 K/2 ±Tk± 7r/a for U >0. The case 
with attractive interaction is shown in Fig. |4ja) with 
the same spread in K as in Fig. [3j In an experiment 
the difference between the main momentum peaks of the 
two species is 2Tk' The case with repulsive interaction is 
shown in Fig. ^h). Alternatively, if the species are held 
by different lasers, the lattice may be turned off adiabat- 
ically in a controlled fashion that keeps the difference in 
the tunneling amplitudes Ja, and hence T^, fixed. 

In this paper we have presented an analytical descrip- 
tion of both the dynamics and the bound states of an 
atom pair of different atoms in a quasi-one-dimensional 
lattice. We suggest several experiments to probe our 
model. The transmission coefficient, Eq. (13), can be 
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FIG. 4: Momentum probability distribution in units of 
{a/2nf for a pair of (a) attractively bound and (b) repul- 
sively bound atoms with the same parameters as in Fig. [3] 
The solid lines show the projected probability of each species 
P{kA,B) in units a/27v. 



a strongly bound state to the structured continuum, and 
the wave number shift Tk of the bound states can be 
measured from the momentum distribution in a time of 
flight experiment. 

We have assumed the two atoms to be in the same 
Bloch band, but they could just as well be in different 
Bloch bands, which also gives rise to different tunnel- 
ing amplitudes J a Jb- In addition, the model can be 
extended to describe two-channel magnetic Feshbach res- 
onances as outlined for identical particles in Refs. [9l [12] 
by replacing Ek and by modifying the Green's function 
in those papers in accordance with our Eqs. ([6|, ([9| and 



probed either by a collision experiment or by RF coupling 
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